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Abstract. The high-order statistics of fluctuations in velocity gradients in the 
crossover range from the inertial to the Kolmogorov and sub-Kolmogorov scales are 
fvq | studied by direct numerical simulations (DNS) of homogeneous isotropic turbulence 

C ■ with vastly improved resolution. The derivative moments for orders < n < 8 are 

t ~— p , represented well as powers of the Reynolds number, Re, in the range 380 < Re < 5725, 

(~~) ' where Re is based on the periodic box length L x . These low-Reynolds-number flows 

\& • give no hint of scaling in the inertial range even when extended self-similarity is 

applied. Yet, the DNS scaling exponents of velocity gradients agree well with those 

deduced, using a recent theory of anomalous scaling, from the scaling exponents of the 

longitudinal structure functions at infinitely high Reynolds numbers. This suggests that 

I~2 \ the asymptotic state of turbulence is attained for the velocity gradients at far lower 

^ ■ Reynolds numbers than those required for the inertial range to appear. We discuss 

K^ , these findings in the light of multifractal formalism. Our numerical studies also resolve 

I— i ' the crossover of the velocity gradient statistics from the Gaussian to non-Gaussian 

behaviour that occurs as the Reynolds number is increased. 
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1. Introduction 

1.1. Motivation and Previous Work 

A deep understanding of the turbulent flow field u(x, t) remains a challenging problem. 
Extensions of the classical theory of turbulence by Kolmogorov [1] consider a multiplicity 
of algebraic scaling exponents for moments of velocity increments 5 r u in the inertial 
range of length scales r, which are spanned roughly between the Kolmogorov dissipation 
scale t]k and the outer scale of turbulence L. The longitudinal increment moments (or 
structure functions) are then given as 



r \ n ( r \ C 



S n (r)E(Ur=^(u(x + r)-u(x)).-j = Ml) ' (1) 

where the scaling exponents ( n depend nonlinearly on the order n but not on the 
Reynolds number i?e, as long as the latter is sufficiently large. The dimensional 
coefficients A n depend at most on large-scale quantities. This nonlinear dependence 
of the algebraic scaling exponents ( n on the moment order n is a manifestation of 
the inertial-range intermittency, which is generally agreed to be an important feature 
of three-dimensional turbulence. Inertial-range intermittency was experimentally first 
quantified by Anselmet et al. [2]. Starting with the work of Kolmogorov [3] and 
Oboukhov [3], numerous phenomenological models have been developed to study and 
describe intermittency (see, for example, Ref. [5]). The most dominant underlying theme 
of these models has been the multifractal formalism [B|. 

There is a similar intermittency in the dissipative scales. This, too, has been 
experimentally characterized since Ref. [TJ, and many models have been developed as 
well (see, again, [5] for a summary). The relation between the two intermittencies has 
been the subject of the so-called refined similarity hypothesis (RSH) put forth in [3]. 
This hypothesis links the statistics of the velocity increments at inertial scales with that 
of the velocity gradients at smaller scales where inertial and viscous ranges match. 

In Refs. O E], it was recognized that dissipation intermittency implies an infinite 
number of dissipative scales, rj. Using this insight, Nelkin [9] worked out the Reynolds 
number dependence of the moments of velocity derivatives. Frisch and Vergassola 
[TO] denoted the range of dissipation scales spanning between r] min and r\ max as the 
intermediate dissipation range. Their geometric picture of the continuum of dissipation 
scales is that each element of the range would possess a local Holder exponent h, 
which characterizes the spatial roughness of subsets of velocity increments in the inertial 
range. Consequently, the minimum and maximum values of the dissipation scale would 
be controlled by the smallest and the largest Holder exponents: r\ min = rj{h min ) and 
Vmax = v(^max)- Later, Chevillard et al. [llj studied the intermediate dissipation 
range within a random cascade model that takes 5 r u as a product of a Gaussian 
random variable and a positive (scale-dependent) transfer variable. They found that 
\n-(Vmax/Vmin) ~ vTmRe. The relation of the intermediate dissipation range to the decay 
of energy spectra was discussed recently in the context of well-resolved shell models [12J. 
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Efforts have also been made to obtain ( n (n) directly from the Navier-Stokes 
equations but the problem has remained a great challenge. We limit ourselves here 
to citing the work of [13] — in part because of the connection to the present work and in 
part because the author kept his considerations close to the dynamical equations. The 
theory has been extended [TU [15] to explore the connection between the viscous and 
inertial range intermittencies. This extension builds on the notion that the fluctuating 
dissipation scale rj is to be considered a field that varies in space and time. A relevant 
feature of the theory is its prediction for the scaling of velocity gradients in terms of the 
exponents ( n . 

Within this overall framework, the present paper accomplishes the following 
goals. First, we perform direct numerical simulations (DNS) of homogeneous isotropic 
turbulence with vastly better spectral resolution than anytime earlier. As was explicitly 
stated by Nelkin [9], such superfine resolutions are required to compute the derivatives 
accurately. Second, we then study the relation between the inertial and dissipative 
regions within the framework of existing theories, namely RSH and theory of [TU [15]. 
After describing the theoretical basis [TJ1[T5] and the details on the numerical simulations 
(Sec. 1), we discuss the analyticity and scaling of the velocity increment moments in 
Sec. 2 and present our findings on the velocity gradient statistics in Sec. 3. We also 
compare in Sec. 3 our results with those of previous work, study in detail the crossover 
of the statistics of the velocity gradients from Gaussian to non-Gaussian regime and 
discuss our results in the light of the multifractal formalism. 

Perhaps the most surprising result of the present work is that, while we find no 
evidence for the inertial range in the DNS data (even when examined through the 
extended self-similarity, or ESS), the measured scaling exponents of velocity gradients 
agree well with those deduced from the longitudinal structure functions at infinitely high 
Reynolds numbers. This suggests that the asymptotic state of turbulence is attained for 
the velocity gradients at far lower Reynolds numbers, well short of those required for 
the inertial range to appear. We conclude with a summary and an outlook in Sec. 4. 

1.2. Numerical Simulations 

The Navier-Stokes equations for an incompressible Newtonian fluid u(x, t) are solved in 
a periodic box of side length L x = 2n. The pseudospectral method is applied with a 2/3 
de-aliasing for the fast Fourier transforms. Advancement in time is done by a second 
order predictor-corrector scheme. The equations are given by 

t? + (u- V)u = -Vp + vAu + f . (2) 

at 

The kinematic pressure field is p(x, t) and v is the kinematic viscosity. We consider 

flows that are sustained by a volume-forcing f (x, t) in a statistically stationary turbulent 

state. This driving is implemented in the Fourier space for the modes with the largest 

wavenumbers kf only, i.e. kj 1 « L x . The kinetic energy is injected at a fixed rate 6j n 
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Table 1. Parameters of the direct numerical simulations. Here, v is the kinematic 



viscosity, £ is the mean energy dissipation rate, R\ = J b / (i£v)u1 ms is the Taylor- 
microscale Reynolds number. We will use the following definition for the large 
scale Reynolds number: Re — u rms L x /v where the box size L x = 2-7T is taken. 
Urms — (u x + Uy + ul) 1 ' 2 is used instead of (#z,u)rms • The spectral resolution is 
indicated by k max rjK where k max = -^ — and N the number of grid points in each 
direction of the cube. 

Run No. N v £ u rms L R\ Re k max r]K 



1 


512 


1/30 


0.1 


0.687 


1.018 


10 


129 


33.56 


2 


1024 


1/75 


0.1 


0.808 


0.920 


24 


380 


33.56 


3 


1024 


1/200 


0.1 


0.854 


0.758 


42 


1074 


15.93 


4 


1024 


1/400 


0.1 


0.892 


0.694 


65 


2243 


9.6 


5 


2048 


1/400 


0.1 


0.882 


0.690 


64 


2218 


19.2 


6 


2048 


1/1000 


0.1 


0.911 


0.659 


107 


5725 


9.6 



into the flow. The volume forcing is established by the expression [T6l [T7] 

f (k, t) = e in ^ ' 5 k , k/ , (3) 

2^k f eK l u l K /> l )\ 

where the wavevector subset K contains kj = (1, 1, 2) and (1,2,2) plus all permutations 
with respect to components and signs. This energy injection mechanism prescribes 
the mean energy dissipation rate; that is, the magnitude of the first moment of the 
energy dissipation rate field, S, is determined by the injection rate, 6j n , having no 
Reynolds number dependence. This can be seen as follows. Given the periodic boundary 
conditions in our system, the turbulent kinetic energy balance, which results from 
rewriting (J2J) in the Fourier space, reads as: 



-V 



^P|u(k,t)| 2 + ^f(k,t)-u*(k,t). (4) 

k k 

The first term on the right hand side of (jl]) is the volume average of the energy dissipation 
field. Additional time averaging in combination with ([3]) results in 

uJ2k 2 (Hk,t)\ 2 ) t = £ = e m = £<f(M)-u*(k,t)) t . ( 5 ) 

k k 

The applied driving thus allows a full control of the Kolmogorov scale i]k = u 3 ^/£ 
in comparison to the grid spacing. In contrast to the usually applied stochastic forcing, 
the integral length scale L, which is defined [18] as 



* r *^, (6) 



2v? x Jo " k 

decreases with increasing Reynolds number. Since the forcing scale in the computations 
is proportional to the box size, we use the box size L x as the relevant scale. The use of 
the integral scale instead of L x does not alter the scaling results significantly Further 
details on the simulation parameters can be found in table [H 
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Figure 1. Resolution requirements in the numerical simulations. The probability 
density function (PDF) of the energy dissipation field £ (x, t) — \(duij dxj + duj / dxi) 2 
is plotted. The dissipation field is given in units of the mean energy dissipation rate 
£. The case with k max r]ji = 1.2 (cyan curve), corresponding roughly to the standard 
resolution in a box of size N = 128, is compared with that of superfine resolution (blue 
curve, see also table [1]). While the cores of both PDFs agree, deviations are manifest 
in the far tails. Both runs are for R\ = 65. Approximately 1.7 x 10 8 data points were 
processed for the analysis in the low-resolution run; the corresponding number for the 
high-resolution run was about 30 times larger. 



As already pointed out, the adequate resolution of the analytic part of structure 
functions turns out to be very demanding. The resolution in the present simulations 
exceed those of conventional simulations by a factor O(10). Consequently, the Reynolds 
numbers attained are modest despite the relatively large computational box, very much 
in the spirit of [19]. In order to stress this point further, we compare the statistics of 
the energy dissipation field £ (x, t) = ^(dui/dxj + duj/dxi) 2 for our resolution and the 
standard case (see figured]). At least for the intended analysis of higher-order gradient 
statistics the proper resolution of the far tails turns out to be necessary. 



1.3. Theoretical Basis 



The theory [TJ1 [15] starts with the exact equations for the n-th order longitudinal 
structure functions S n {r) which can be directly derived from the equations of motion 
for the turbulent fluid PS] [20]. For homogeneous, isotropic and statistically stationary 
turbulence in three dimensions these equations take the form 

2(2n- 1) 



or r 



-G 2n -2,2(r) + (2n - l)S r a(8 r u) 



2n-2 



(7) 



Here, G2 n -2,2 is the mixed term containing longitudinal increments of order 2n — 2 and 
transverse increments of order 2. Equation ([7]) is not closed because the last term on the 
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Figure 2. The matching of the inertial and dissipative ranges, (a) In the standard 
approach, the same Kolomogorov scale matches the singular (~ r^ n ) and analytic parts 
(~ r n ) for all orders n of the increment moments, (b) In the theory of [HI [T5], for 
each order of the moment S n (r), the singular and analytic parts match at an order- 
dependent scale r\ n . Following Q12j). r\ n < r\ m when n > m. 



right hand side is unknown. For small increment scales around the Kolmogorov length, 
it follows that the expression for the Lagrangian acceleration of fluid particles is given 

by 

(<V0 3 



Srjd 



(8) 



recalling that the characteristic time is of the order v/{8 v u) 2 . In this equation, note 
that 7] is a field and that increments are therefore taken across variable distances. The 
unknown term has the form 

1- 



5 r a(S r u] 



2n-2 



V 



-(5 v u) 3 (5 r u] 



2n-2 



(9) 



This correlation involves two scales — a locally varying dissipation scale field rj in the 
acceleration increment and r in the velocity increment moment — and is therefore hard 
to manipulate. However, in the limit r — > n one can make some progress. In this limit, 
we set rj = r]2n, where r]2 n is the order- dependent matching distance between the analytic 
and singular parts of S l 2 n (r) (see figure [2]). The other three terms of (j7j) are found to 
be of the same magnitude (2TJ, [22] and are of the order S^nj^in- One thus obtains, for 
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r — >■ 77, the result 

S2n{l]2n) ^ SWfll^n) /-, „n 

??2n f 

Now, the velocity increments have the property that, at the large scale L, their 
distribution is Gaussian. It then follows from (1) that Sznivin) = (2n— l)\\a]^ (r]2 n / L)^ 2n , 

where <jl = y (o~lu) 2 - lil Putting r = r] 2n , one obtains from (5) 



/7/2n\ Can ~ fat + 1 ~ 1 _ (2n)!! a L L _ a L L 

V L / ~ (2n-l)!! u ™ v ' l j 

With the large scale Reynolds number i?e = a^Ljv one gets 

772n « Iife^-^1- 1 . (12) 

For the Kolmogorov scaling, ( n = ~ and (TT2l) yields r] 2n = LRe~ 3 ^ 4 = tjk for all orders 
n, as consistency would require. 

To make further progress, the functional dependence of C,2n has to be given explicitly. 
The theory of [T3j provides a convenient functional form 

2(l + 3/?)n 
C2n "3(l + 2/5n)' [i6) 

which, with the free fitting parameter (3 set to 0.05, agrees with available experimental 
data in high- Reynolds-number flows (for order 10-15). We find that this relation 
agrees, up to order 10 — 15, with available measurements as well as with popular 
parametrizations of ( n , e.g. with the She-Leveque model and the p-model [231 121] ■ 
The scaling behavior of the spatial derivative in the analytic range of the displacement 
r can be calculated, in the limit r — > 77, as 



du 



dx 



5 r] u 



V 



(S v u 



2n 



o^Re n 4n=Re pn , (14) 



v" 

where we have used that the "dynamic" Reynolds number at the local dissipation scale 
is unity, i.e. Re v = r]8 v u/v ~ 1. The use of (TP2"|) yields 

Pn = n + -^ . (15) 

S,2n — C,2n+1 — 1 

Since (3 = 1 [25]) relation (|T5"|) gives (2 = (2 — 2pi)/(2 — pi). For the Kolmogorov 
value of (2 = 2/3, we obtain p\ = 1/2. The anomaly that may exist in the first-order 
exponent p\ for velocity gradients is related to the second-order inertial exponent £2- 
For moments of the dissipation rate, one can write 



where 



d n = n + - ^ -. (17) 



| For example, S4,(L) — 3!!cr^(L/L) l > 4 = ia\ and S%{L) = a\. It follows that, at the appropriate large 
scale L, the flatness factor S^L)/ Sf (£) = 3. 
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Thus, the second-order quantities for the dissipation rate are expressed in terms of the 
eighth-order quantities involving velocity increments. In general, to accurately evaluate 
£ n , one has to resolve the analytic range within which 5 < 4 n (r) oc r 4n . This difficulty for 
large n is one of the main considerations of the theory. 
Finally, we note that 

di = (18) 

because £ = €i n holds in the DNS (see equation (jSJ)). From relation (|T7j) one immediately 
has 

Cs = 2C 4 - 1 • (19) 

Measurements are, in fact, in good conformity with this equation. 

l.Jf.. Velocity Gradients from Refined Similarity Hypothesis 

The refined similarity hypothesis (RSH) [3] imposes a different constraint between 
moment orders of the energy dissipation and structure functions. This results in a 
different Reynolds number dependence compared to (IT7|) . When taking £ ~ (5 Vk u) 3 /vk 



(20) 



(21) 

The comparison with the present data will be made later in this text (see table [2]). Here, 
we briefly mention that the intermittency exponent // in the scaling £ (x + r)£(x) ~ r~ M 
is 0.25 from RSH while the application of the theory of Refs. [HJ [15] gives /i ~ 0.2. 
Both are within the accepted range of 0.25 ± 0.05. 

1.5. Velocity Gradients from the Multifractal Formalism 

Nelkin [9] predicted a Reynolds number dependence for derivative moments based on 
the multifractal formalism (MF). The derivation of his expressions relied on the steepest 
descent approximation of resulting integrals and input from measurements [TJ. His result 
(which is also outlined in detail in [6], pp. 157-158) is that 

p n =p(n)-n, (22) 

in terms of the notation used above. Here p(n) is a unique solution that follows from 
the intersection of the concave curve ( p with the straight line 2n — p for a given n. With 
(TT3|) the intersection (p > 0) obeys the relation 

p( w ) M _ 4 + 3 ^ ~ 2n ) + --^72/fa + (3/?(2n - 1) - 4)2 . (23) 

The resulting values of p n are also listed in table [2j 



the relation 












£ n ~ 


S3n(VK) 

Vk 


-Vk" 


n 


follows. With 


Vk = 


: LRe- 3 / 4 

_ 3(n-C3„; 

Re 4 


we get 
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Figure 3. Longitudinal structure functions of the turbulent velocity field, (a) Second- 
order longitudinal structure function Siiy^jv 2 ^ over r/?7x for the three different Runs 
indicated in the legend, (b) Tenth-order longitudinal structure functions Sio(r)/v 1 ^ 
over r/r)K for the same data. The inset shows that the curves do not collapse well. 
Here, % = (V5) 1 / 4 is the Kolmogorov velocity, (c) The tenth-order structure functions 
for R\ = 24, 42, 65 collapse when r is rescaled by the dissipation scale tjiq defined by 
(fT2| and the amplitudes by the velocity scale v\q (see equation (|25|)). The inset shows 
the same level of expansion as in (b). (d) The figure shows the exponent p n as a 
function of order n given by equations (|f3[) and (If 5[) . 



2. Moments of Velocity Increments 

2.1. Rescaling of Higher Order Moments and Test of Analyticity 

Figures Ek and [3b show the structure functions S2 n (r)/v]^ plotted against r/r/x for three 
Reynolds numbers for order 2 (n — 1) and 10 (n = 5). Here, vk is the Kolmogorov 
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10 




Figure 4. Test of analyticity for the longitudinal structure functions S-2 n {r) by 
compensated plots. Left: Orders 2n = 2,4,6,8,10 are shown for Run 5. Horizontal 



lines indicate the exact analytical form corresponding to S2n(r) 



Right: Sio(r) 



is shown for Runs 2, 3 and 5. The red vertical arrows indicate the corresponding 
Kolmogorov scales t\k ■ 



velocity corresponding to the definition i]kVk/v = 1- Two features of the graphs are 
worth noting. First, no inertial range can be seen at these low Reynolds numbers. 
Second, all data possess the analytic parts thus confirming that the resolution used 
is adequate. It is important to recognize in figure H] that with increase of both the 
Reynolds number and moment order the width of the analytic range decreases. To 
make this point more explicit, we plot in figure H] the compensated structure functions 



Sin/r 



2n 



n 



5) for a fixed Reynolds number R\ = 64 (left) and the normalized 



moment Sio/r 10 for different Reynolds numbers (right). In the analytic range we expect 

S 2n (r) 



,-2n 



fdu\ 

\ dx ) 



2d 



const . 



(24) 



i.e., independent of the increment distance r. For R\ = 64, this range is well-defined 
for the moments n < 3 and just survives for n — 5. This means that, even for the 
present super-resolution, the representation of moments of velocity derivatives in terms 
of the low-order finite differences may be problematic for moment orders higher than 5 
(although we do present the data for n up to 8). 

In the analytic range, all curves can be expected to collapse when normalized by 
the appropriate length and velocity scales. The traditional scales are the Kolmogorov 
length and velocity scales rjK and vk, respectively. This scaling works well for low order 
moments (say 2), as seen in figure [3k. The same normalization is not adequate for high 
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orders such as 10, as can be seen in figure [3b. However, all curves do collapse when the 
length scale r] 2n (see equation (fl2l) ) and the corresponding velocity scale 

V 2 n = — (25) 

1l2n 

are used instead of r] K and vk (see figure [3b)l§| 

How can we understand this collapse? Introducing the "dynamic" Reynolds number 
Re r = r5 r u/u, and recalling that the only relevant parameter in the inertial range — 
including the interval just above the dissipation range (r > 77) — is the energy flux, we 
can conclude that the dynamics of fluctuations at the scales r ps r\ are independent 
of the width of the inertial range (L,rj). While one cannot, in principle, rule out the 
existence of other fluxes corresponding to some yet unknown integrals of motion related 
to structure functions of higher order than 2 [26J , this possibility will not influence the 
essence of the argument. If so, for r < 77, the properly normalized moment of a given 
order n must be independent of Re. 

2.2. Test of Extended Self- Similarity 

Since no inertial-range scaling can be observed for such low Reynolds numbers in the 
standard double logarithmic plot, we tested the method of extended self-similarity (ESS) 
by Benzi et al. [27J. As is well known, the main point of ESS is that even though there 
is no discernible scaling in the standard plot when the Reynolds numbers are low, one 
can detect a sizeable range of scaling when arbitrary moment orders are compared with 
S 3 : 

S n (r) ~ (-S 3 (r)) c » • (26) 

The sensitive and scale-dependent measure for possible anomalous scaling is then the 
local slope 

dln(g n (r)) 
Xn(r) = , , „ • 27 

dln(-5 3 (r)) 

The results for our data are shown in figure [5] for the fourth and tenth order. In both 
figures we cannot detect anomalous scaling in the ESS framework, which would result in 
local slopes of Xi( r ) ~ 1-28 and Xio( r ) ~ 2.59 [23J[T3]. We conclude that only a further 
increase of the Reynolds number will shift the local slope toward the asymptotic values. 
Thus, even the "backdoor" of ESS is not opened for our small Reynolds numbers. It 
confirms our statement made before that the velocity field statistics does not have any 
asymptotic scaling in the inertial range. The inset reveals an interesting feature. It can 
be observed that a faster relaxation towards anomalous scaling occurs with increasing 
order. At this point, we can only speculate about the reason for this feature. On one 
hand, the differences between the viscous scaling, Xn( r ) = n/3, and the inertial scaling, 

§ We note here parenthetically that S10 for the lowest Reynolds number (R\ = 10) does not collapse 
on the common curve. For this case, however, the velocity gradient statistics are Gaussian, in contrast 
to the other three cases. This issue will be discussed in greater detail in Sec. 3.2. 
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Figure 5. Extended self-similarity analysis of longitudinal structure functions, (a) 
Fourth-order structure function ^(r) versus — Ss(r). Structure functions are defined 
by equation (QJ, with — S^r) = (4/5)£r. (b) Tenth-order structure function Sio{r) 
versus — ^(r). The insets in both figures show the corresponding local slope Xn(r) as 
given by relation (f27|) . Data of Runs 2,3 and 5 are indicated by different colours, as 
noted in the legend. 



x{r) = ( n , become larger with increasing order and thus better visible in the local 
slope. On the other hand, an order-dependence of the dissipation scale, n n , might cause 
a slight increase of the inertial range and the crossover between inertial and viscous 
ranges, respectively. 



3. Velocity Gradient Statistics 

3.1. Results and Comparison with Refined Similarity Hypothesis 

If indeed the properties of fluctuations from the interval r < n depend only upon the local 
magnitude of the energy flux S and not upon the width of the inertial range (L, rj), one 
can hope to obtain the asymptotic values of exponents p n (Re) -^ p n (oo) in reasonably 
low-Reynolds-number flows, provided that a small (even very small) constant-flux range 
exists for scales r > n. 

In figure El we plot \du/dx\ and £ n as functions of Re. Statistical convergence 
of both gradient quantities is satisfactory for gradient moments of at least up to the 
7-th order (see figure [6b) and dissipation moments of at least up to the 4-th order (see 
figure (6U). As can be seen from table EJ the data agree well with theoretical predictions 
[14"l [15] . For example, we get p\ = 0.455 compared to the theoretical value of 0.465. 
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Figure 6. The Reynolds number dependence of the moments of velocity derivative 
and energy dissipation, (a) Moments of the absolute value of the longitudinal velocity 
derivative du/dx for orders 1, 3, 4, 5 and 7 as functions of the box Reynolds number 
Re. Only the filled data points were included in the least square fit, but Run 6 
demonstrates that the scaling continues for higher Reynolds number. The case with 
the lowest Reynolds number has Gaussian statistics (as will be described below) and 
is hence not turbulent in the traditional sense, (b) Moments of order 2, 3 and 4 of 
the energy dissipation field £. Again, only the filled data points were used for the 
fit to evaluate the exponent, (c) Statistical convergence test for the seventh order 
longitudinal velocity derivative moment at the highest Reynolds number considered 
here, R\ = 65 (Run 4). The data set contained 15 samples of the turbulent field which 
results in 1.6 x 10 10 data points, (d) Statistical convergence test of the fourth order 
moment of the energy dissipation field. Again, data from Run 4 are used. 



The second order exponent (2 is given by 
2-2p! 



<2 



0.706 



(2* 



2-pi 

which is very close to the experimental value of (2 — 0.71 [271 [23]. We stress that this 
result was obtained in flows with 24 < R\ < 65, none of which has any inertial range. 
The moments of the dissipation rate also agree with theoretical predictions, as shown 
in table [3 
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Table 2. Comparison of scaling exponents for different velocity gradient moments. p n 
for n = 1,3,5,7 (see (|15l0 and d„ (see (fT?} ) for n = 1,2,3,4 are listed. Results from 
the present DNS are compared with those from the theory [21 |T5] after inserting (fT3")) 
into (fTS"]) and (|17|) . Comparisons with the refined similarity hypothesis (RSH) and the 
multifractal formalism (MF) are also provided. The error bars for orders 3 and 4 of 
the DNS data have been determined from lower and upper envelopes to the tails of 
£ n p(£), p(£) being the probability density function of the energy dissipation field. The 
same holds for orders 5 and 7 in case of p(\du/dx\). The range of Reynolds numbers 
for all fits was 380 to 2243. 





Theory [14] 


DNS 


RSH 


MF 


Pi 


0.465 


0.455 


- 


0.474 


p?, 


1.548 


1.478 


- 


1.573 


Pi 


2.157 


2.051 


- 


2.188 


P5 


2.806 


2.664 ± 0.137 


- 


2.841 


P7 


4.203 


3.992 ± 0.653 


- 


4.241 


<*i 


0.000 


0.000 


0.000 


- 


d 2 


0.157 


0.152 


0.173 


- 


d 3 


0.489 


0.476 ± 0.009 


0.465 


- 


di 


0.944 


0.978 ± 0.034 


0.844 


- 



In the theory outlined earlier [H], the second, third and fourth moments of £ are 
related to the structure functions S 8 (r] 8 ), 512(7712) an d Si e (r]i 6 ), respectively. These 
structure functions probe very intense, low-probability velocity fluctuations and, as a 
result, huge data sets are needed to accurately evaluate their characteristics. We have 
seen that dissipation moments of order 5 are barely resolved in the present simulations, 
while the data for the sixth-order moment S 6 , corresponding to S24 have not converged 
well. In addition, we stress that the statistical convergence is not sufficient for the 
accurate determination of structure functions: the simulations must resolve accurately 
at least a fraction of the analytic range r < r\± n . 

Table 2 lists comparisons between theoretical considerations and the DNS data. 
In case of the dissipation field, we inserted ( TT3l) into ( 1211 ; for the velocity derivative 
scaling exponents from multifractal formalism, we used the relation (J23I) which was 
inserted into ( J22l . For the gradient exponents, the DNS data are somewhat smaller 
than both the theory [T4"l [15] and the multifractal theory. In the case of the dissipation 
exponents, the DNS results are closer to the theory of Refs. [i~4l [j~5] and somewhat 
larger than those of RSH. Note that the results of both theories depend on the inertial 
scaling exponents obtained from measurements (or at least expressions tuned to agree 
with measurements). Considering that there are issues of resolution in measurement 
(although not in the same sense as in simulations), these departures may suggest that 
inertial exponents may need slight revision. Of course, there might be other reasons for 
these differences. 
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Figure 7. Statistics of the velocity gradient du/dx. (a) Plots of the probability density 
functions (PDF) for Runs 1 to 4 as indicated in the legend. The gradient is normalized 
by its root-mean-square value. For comparison the Gaussian distribution is added to 
the figure, (b) Corresponding skewness of the PDFs (blue symbols), (c) Corresponding 
flatness of the PDFs (blue symbols) . In order to highlight the transition, further data 
points have been added to the data of table Q] in figures (b) and (c) . These additional 
numerical simulations have been conducted at a spectral resolution of N — 256 for 
different kinematic viscosities. Forcing scheme and energy injection rate are the same 
as before. The additional data values are plotted as cyan filled squares. 



3.2. Transition from Gaussian to non- Gaussian Velocity Gradient Statistics 

As mentioned earlier, while computing scaling exponents, Run 1 was excluded from the 
least-square fits. The reason is the qualitatively different nature of the velocity gradient 
statistics at this lowest Reynolds number (R\ = 10). This will be discussed now. Figure 
[7] illustrates the crossover from Gaussian statistics at very low Reynolds numbers to 
an increasingly non-Gaussian behaviour for moderate Reynolds numbers. In order to 
highlight this transition, we generated additional DNS data at intermediate Reynolds 
numbers. The Gaussian values for the third and fourth order normalized derivative 
moment are indicated by the red solid lines in panels (b) and (c) of figure 7. The 
derivative moments are defined as 



M n {du/dx) 



(du/dx) n 

W2 ' 



(29) 



(du/dx) 2 

For the lowest Reynolds numbers we detected a regime which is a complex time- 
dependent flow rather than a turbulent one. In this regime, the flow can be basically 
described by the small number of driven modes that form a low-dimensional nonlinear 
dynamical system. All other degrees of freedom are strongly damped and slaved to the 
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Figure 8. Morphological manifestation of the crossover from Gaussian to non- 
Gaussian velocity gradient statistics. Isolevel plots of snapshots of \du/dx\ for Run 1 
at R\ = 10 (top) and for Run 4 at R\ = 65 (bottom) are shown. Both level sets were 
taken at the corresponding values of 2 x (du/dx) rms . We observe a significant increase 
of the spatial intermittency for the higher Reynolds number data set. 
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driven modes. The sign of the derivative skewness M 3 (du/dx) there became sensitively 
dependent on the particular modes that were driven. Since the non-Gaussian behavior 
is related to the acceleration- velocity term in equation (7), the experience suggests that 
it is difficult to pin down the behavior of this term near zero. The magnitudes of 
M3 and M4 varied, respectively, around zero and three (even for very long-time runs). 
We verified this by choosing different wave vectors for driving the flow while leaving 
all other simulation parameters the same, including the energy injection rate e« n (see 
equation (J5])). At larger Reynolds numbers, R\ > 10 or 15, the derivative moments 
became insensitive to the particular driving. Therefore, our studies suggest that the 
transition to non-Gaussian statistics is smooth with respect to the Reynolds number 
and that intermittency continuously grows with the growing number of excited modes. 
Figure El illustrates the morphological changes of the spatial distribution of the 
velocity gradient that are connected with the change of the statistical properties. The 
increasing intermittency of the velocity gradients is accompanied by an increasing 
fragmentation of the isolevel sets. 

4. Summary and Discussion 

Turbulence in a three-dimensional periodic box, generated by the Navier-Stokes 
equations driven by a large-scale forcing, was investigated in the Taylor microscale 
Reynolds numbers range 10 < R\ < 107. The simulations were made with superfine 
resolution in order to resolve the analytic part of structure functions at least up to order 
16. No inertial ranges characterized by the velocity structure functions S n (r) oc r^ n 
were detected — not even with the method of extended self-similarity In the range 
24 < R\ < 107, strong intermittency of the spatial derivatives was detected and 
their moments were accurately described by the scaling relations \du/dx\ n ex Re Pn 
and £ n oc Re dn , respectively. The exponents p n and d n were found to be in essential 
agreement with the theoretical work for high Reynolds number [HI [T5] . Based on the 
well-resolved numerical results, we are thus able to relate the dissipation-range exponents 
of the velocity gradients to the inertial range scaling exponents of the velocity field for 
very high Reynolds numbers. For instance, the DNS result p\ ~ 0.455 gives £2 = 0.706 
corresponding to the asymptotic high-Reynolds number energy spectrum E(k) oc &;~ L706 . 

The competing theoretical predictions are from RSH (for dissipation) and the 
multifractal theory (for the gradients). They, too, agree with the DNS results, although 
somewhat less successfully overall. It appears that the theory of Refs. [T4"l [j~5] seems to 
have an edge. 

This last conclusion would have been more conclusive if we had been able to obtain 
reliable dissipation range statistics accurately for moments of order 5 or 6. Despite the 
huge data sets of ~ 10 10 points, we were unable to obtain reliable data for dissipation 
field moments with n > 4, certainly for n no larger than 5. The theory suggests the 
reason for this problem: moments of the dissipation rate S n are related to high-order 
structure functions S in {r}^ n ) probing very low probability fluctuations evaluated on the 
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corresponding dissipation scales. 

We also conducted simulations at various Reynolds numbers to study the smooth 
transition of the velocity gradient statistics from the Gaussian behaviour to the non- 
Gaussian behaviour. 

Our work suggests that the scaling exponents of the moments of velocity derivatives 
observed in the relatively low- Reynolds-number turbulent flow, lacking even traces of the 
inertial range, can be expressed in terms of the inertial-range exponents corresponding 
to the asymptotic case (Re — ► oo). We stress that the scaling exponents p n and d n of 
the moments of velocity derivatives and dissipation rate reach asymptotic values that 
are independent of large-scale Reynolds number, even at low values of the Reynolds 
numbers. Thus, the dynamics of velocity fluctuations at the scales r m i] are asymptotic 
even in relatively low-Reynolds-number flows. This could mean that the magnitudes of 
inertial-range exponents (Re — > oo) are prescribed by the matching conditions on the 
ultra-violet cut-offs formed in the low-Reynolds-number regimes. The relation of the 
observed behavior to fluctuations of the dissipation scale will be discussed elsewhere. 
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